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In this paper we introduce an equation of state (EOS) of quark matter within the framework of 
Dyson-Schwinger equations (DSEs) to study the structure of compact stars. The smooth crossover 
from hadronic matter to quark matter in the hybrid star is studied. We compare different strategies 
to obtain crossover EOSs and find a new way to construct two-solar-mass hybrid stars with even a 
relatively soft quark EOS, while earlier works show that the quark EOS should be stiff enough to 
support a massive hybrid star. 
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I. INTRODUCTION 

The equation of state (EOS) of quark matter is crucial 
for the study of quark stars and compact stars. Once the 
EOS is obtained, the radius and gravitational mass of a 
compact star with a given central mass density can be 
calculated by solving the Tolman-Oppenheimer-Volkoff 
(TOV) equations. By changing the central density, the 
mass-radius relationship curve that is used to be com¬ 
pared with the astronomical observations can be drawn. 
Thus, the astronomical observations can be used to con¬ 
strain the parameter sets of the EOSs and even rule out 
some EOSs. However, on the one hand, the perturbative 
QCD method fails in the low-energy region, and it is very 
difficult to obtain a reliable EOS from the first principles 
of QCD. Some approximate methods that incorporate ba¬ 
sic features of QCD are adopted to calculate the EOS of 
quark matter. For example, the Massachusetts Institute 
of Technology (MIT) bag model and the Nambu-Jona- 
Lasinio (NJL) model. But all these models have their 
own weaknesses. The MIT bag model violates chiral sym¬ 
metry even in the limit of a massless quark. So, there are 
also some authors developed advanced bag models such 
as the so called chiral MIT bag model [T]. This model 
restores the chiral symmetry by introducing the coupling 
of pion and quark. The NJL model assumes that the 
interaction between quarks is point-like, so this model is 
not renormalizable, and it cannot incorporate quark con¬ 
finement [2]. On the other hand, the observational max¬ 
imum mass of compact stars is increasing. For example, 
the P5'RJ0348 + 0432 was reported to have 2.01 ± 0.04 
solar mass in 2013 [3]. Only very stiffen EOS can sup¬ 
port such a large maximum mass. Nonetheless, the EOS 
is often softened when hyperon mixing, deconfined quark 
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matter or kaon condensate is taken into consideration in 
the core of a compact star. So, F. Ozel even argued that 
the condensates and unconfined quarks may not exist in 
the centers of neutron stars |4]. 

To solve these problems, people tried to find stiffen 
EOSs that can support massive quark stars and neutron 
stars. For instance, in a pure hadronic framework, the in¬ 
troduction of a “universal 3-body force” acting on all the 
baryons or an effective Lagrangian includes quartic terms 
in the meson fields can make the hyperonmixed neutron 
stars become possible I3E]. In MIT bag model, some au¬ 
thors also introduced a density dependent bag parameter 
to get stiff EOSs [7]. Besides, the smooth crossover from 
hadronic matter to quark matter is used to construct hy¬ 
brid stars with high mass recently. Certainly, the order 
of hadron-quark phase transition at zero temperature is 
still an open problem and the first-order phase transition 
(Maxwell or Gibbs construction) is widely adopted even 
in some recent works BM, and the exist of the so-called 
“mass twins” in the mass-radius relationship for compact 
stars is presented as an potential evidence of the first or¬ 
der phase transition m- However, we know that the ra¬ 
dius of a compact star is very difficult to be determined 
exactly, so it is not easy to find “mass twins”, and there 
are also some arguments that the “mass twins” also exist 
in the smooth phase transition case m- Besides, lat¬ 
tice QCD (LQCD) studies showed that the transition line 
for low net baryon is a crossover [12HIS], but the results 
for low temperature and large baryon density are still 
model dependent. Many people suppose the existence of 
a critical end point (CEP), but others do not. Actually, 
on the theoretical side there is still an ambiguity, not 
only for the location of CEP but also for the existence of 
CEP [16] . For example the authors of Ref. m argue that 
if the vector interaction is strong enough, the transition 
is a crossover in the whole phase diagram. Furthermore, 
treating the point-like hadron as an independent degree 
of freedom in the Gibbs condition is not fully justified 
in the transition region because all hadrons are extended 
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objects composed of quarks and gluons. The study on 
smooth phase transition in hybrid stars is thus also nec¬ 
essary. The picture of a gradual onset of quark degrees of 
freedom in dense matter associated with the percolation 
of finite size hadrons has been discussed in some semi¬ 
nal works such as Refs. nil US]. Recently, the smooth 
crossover from hydronic matter to quark matter is used 
to construct massive hybrid stars [2QU22] . In these works, 
based on the percolation picture, the authors utilized dif¬ 
ferent interpolation functions to do some phenomenolog¬ 
ical studies on the hybrid stars compatible with 2 solar 
mass. Obviously, the strategies of interpolation are dif¬ 
ferent, but in each paper, the authors just concentrate 
on a certain interpolation strategy. 

The aim of this paper is to find a proper EOS of quark 
matter to calculate the structure of compact stars and 
discuss the possibility of two solar mass hybrid star based 
on the smooth phase transition between the hadronic 
matter and quark matter. So, we perform a model study 
on the hybrid star from the point of view of the smooth 
crossover. For the quark phase, we adopted the EOS 
with strangeness based on the Dyson-Swinger Equations 
(DSEs), because the DSEs can simultaneously address 
both confinement and dynamical chiral symmetry break¬ 
ing [23l |24] , and it has been applied to hadron physics 
successfully [251130] . We discuss different interpolation 
strategies to construct hybrid EOSs and compare their 
influence on the final results, that is to say, the maxi¬ 
mum mass of the hybrid stars. We find a new way to 
construct an EOS that is sufficiently stiff to include sta¬ 
ble hybrid stars compatible with two solar mass and our 
conclusion is an excellent complement to the conclusion 
in the Ref. [20] . Because it is usually considered that the 
EOS of quark matter should be stiff enough to support 
a massive hybrid star compatible with two solar mass 
before. 

This paper is organized as follows. In Sec. II, we intro¬ 
duce the EOS of quark matter. In Sec. Ill, several ways 
are tried out to construct hybrid EOS from the point of 
view of the smooth crossover phase transition. A new 
way to construct massive hybrid stars even with a soft 
EOS of quark matter is proposed. Finally, a brief discus¬ 
sion and conclusion is presented in Sec. V. 


field-strength renormalization constant, is the mass 
renormalization constant where m is the current quark 
mass, g is the coupling constant. Zip is the quark-gluon- 
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vertex renormalization constant, / ( 2 ^ ^ symbol 

that represents a Poincare invariant regularization of the 
four-dimensional Euclidean integral, is the GellMann 
matrices, r^(p, q) is the quark gluon vertex, and D^y{p — 
q) is the dressed gluon propagator. The general form 
of the quark propagator at zero temperature and zero 
chemical potential reads 

S{p)~^ = ii)A{p^) + ( 2 ) 

where A{p^) and B{p^) are scalar functions of . The 
renormalization condition is 

^ (c") = r (3) 


B (C^) = m, (4) 

at sufficiently large space-like [32l [33]. We choose 
the renormalization point to be C = 19 GeV, as used 
in Refs. [32l|33]. Another popular choice is ^ = 2 GeV in 
recent study, e.g. Ref. [34]. Actually, DSEs is a renormal- 
izable theory, so the choices of C do not affect the physical 
results.. Eq. (1) is the exact result from the first prin¬ 
ciples of QCD, but we can not solve it directly unless 
concrete truncations are performed. Rainbow truncation 
is used in this work, which means a bare vertex is adopted 


ya 

= -jiv, 


( 5 ) 


and the Qin-Chang gluon propagator model [35| is spec- 
ified by a choice for the effective interaction in Landau 
gauge. Finally, one can obtain the following equations 
for the two dressing functions A{p^) and B{p^)^ 




g(tf) A(g^) 


X 


p-g + 2 


{k-p){k • ^)\ 

k^ ) ’ 
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II. THE EOS OF QUARK MATTER AND THE 
DSES 


B (p^) = Z4^m + 4 


L 


g {k^) B (g^) 

k2 g2^2 (^q2^ _|_ ^2 (^q2^ 


( 7 ) 


In this section, we briefly introduce our EOS of quark 
matter based on the DSEs approach. According to Ref. 
m, we start from the zero temperature and zero quark 
chemical potential version of the quark propagator DSE 
that reads 

Sip)~^ = Z2{ii) + z^m) 

+ g^ZiF Y 7 ^S'(g)r“(p, q)D^„{p - q), ( 1 ) 

where S{p) is the dressed quark propagator, Z 2 is the 


The coupled Eqs. (6) and (7) can then be solved by 
direct iteration. The DSE solution for the dressed quark 
propagator can be well fitted with 3 pairs of complex 
conjugate poles with the representation [36| 


s{p) = E 


n=l 


ipl -h ipS -h m* 


(8) 


We fit u, d and s quarks respectively, where the following 
requirement, that the quark propagator in the region of 
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ultraviolet should tend to the free quark propagator, is 
employed 
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i^k + zl) = 1. (9) 

k=l 

In this paper, we choose a new parameter set that can 
not only fit the propagator well but also lead to a stiff 
EOS. The corresponding parameters for u, d quarks are 

mi = 285 + 121i MeV, = 0.308 + 0.618i, 
m 2 = -1120 + 59i MeV, ^2 = 0.111 - 0.193i, (10) 

ms = 1214 + 429i MeV, zs = 0.081, 

and for s quark. 


mi = 509 + 236i MeV, zi = 0.327 + 0.449i, 
m 2 =-1166 + 6162 MeV, ^2 = 0.101 + 0 . 001782 , 
ms = 1572 + 6602 MeV, ^s = 0.072 + 0.0172, 

(11) 

If we compare these parameters with that in the Ref. 
[31], we will find that the smaller real parts of mi and 
zi make the EOS stiffer. Then, the propagator at zero 
temperature and zero chemical potential case can be gen¬ 
eralized to nonzero temperature and chemical potential 
case by the following replacement [37] 

COfi = (2n + l) 7 rT + i/i, ( 12 ) 


which is widely used within rainbow truncation of DSEs 
in thermal and dense QCD, and Cbn is the Matsubara 
frequencies. Hence, the quark propagator at nonzero T 
and p is 


S 




Zk 


+ r- 


^ V ¥ + 174‘I’n + 'mk ^ + i 74 ‘l’n + f^k/ ' 

(13) 

Then, the well-known formula for the quark number den¬ 
sity is 


= -N,NfT 


+ 00 » 

n= — oc 


(P p 



(14) 

After a series calculations and Taking the limit T ^ 0, 
the final result of the quark number density in zero tem¬ 
perature can be obtained: 


p{p,T = 0) 

= E ^ (m - Pk) (1^" - ^ - Cfe) ^ 

(15) 

where = |Re(m/c)| and Ck^dk are defined by m^ = 
Ck + dki. The quark number density dependence on p 
at T = 0 for s quark is displayed in Eig. 1, and from 
this figure we can see that the quark number density of s 
quark turns to nonzero at pc = 520 MeV, which is smaller 
than the Ref. m because of our new parameter set. So, 
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FIG. 1: The quark number density of s quark 



this EOS is more appropriate for the quark matter with 
strangeness. Einally, we take the chemical equilibrium 
and electric charge neutrality conditions into considera¬ 
tion to constrain different quark chemical potentials. The 
conditions read 


Pd — pu T pe^ 

Ps pu T Pei (^7) 

2 11 

-^Pu - -^Pd - -^Ps - Pe = 0 . (18) 

Then there is only one independent chemical potential 
due to these constraint. We choose pu in this work. Eor 
definite chemical potential of quark, the EOS of QCD at 
T = 0 reads [38l 

P{p) = P{p = 0)+ [ p{p')dp', (19) 

Jo 

and the relation between the energy density and the pres¬ 
sure of the corresponding system is gniiiT] 

e = —P + ^^piPi^ (20) 

where pi and pi represents the chemical potential and 
the particle number density for each component in the 
system. 


III. THE SMOOTH CROSSOVER AND THE 
STRUCTURE OF HYBRID STARS 


As usually, one can obtain the structure of a bare quark 
star by directly integrating the Tolman-Oppenheimer- 
Volkoff (TOV) equations with the equation of state: 


dP{r) 

dr 


G{e + P) (M + dTrr^P) 


r(r - 2GM) 
dM{r 


= 47rr 5. 


( 21 ) 


dr 
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However, in the Eq. (19), we can see that the P{ii = 0) 
is only a /r-independent constant. It is the pressure of the 
vacuum. Although a method to calculate the P(/i = 0) 
self-consistently is given in the Ref. m, there is actu¬ 
ally no model independent way to calculate the reliably 
P{lJi = 0) from the first principles of the QCD. In this pa¬ 
per, analogous to the MIT bag model, we reconsider the 
effect of this term and think that there exists negative 
pressure at zero chemical potential in the vacuum which 
manifests the confinement of QCD. Namely, we identify 
P(/i = 0) with -B, where B is the vacuum bag constant. 
In this paper we take it as a phenomenological parame¬ 
ter only. Our EOS for nuclear matter is the APR EOS 
with AI8 P 6v P UIX"" interaction in Ref. [42]. This 
EOS is based on Argonne vI8 two-body potential and 
the Urbana IX three-body interaction, includes charge 
neutrality and Beta-equilibrium. The 6v indicates the 
inclusion of relativistic corrections. Eor simplicity, the 
APR model includes only nucleonic degrees of freedom, 
and does not take into account hyperons, whose interac¬ 
tions with nucleons and among themselves are not well 
determined. Then, we choose B = (110 MeV)^ to en¬ 
sure that the energy density of the quark matter should 
always be higher than the hadronic matter in the low 
density region. But with this value of B, the maximum 
mass of a pure quark star is just about 1.65 solar mass. 
If we choose a smaller B, the EOS will become stiffer 
and the maximum mass will become larger, but such a 
small B is unreasonable. So, we just try to find a way 
to construct massive hybrid stars in the point of view of 
smooth crossover phase transition with such a soft EOS 
for quark matter. In order to facilitate analysis and com¬ 
pare with the first order phase transition, the pressure 
as a function of the baryon chemical potential of both 
the quark matter and the hadronic matter is showed in 
the Eig. 2. Erom Eig. 2 we can see that the intersec¬ 
tion in the P — jiB plane is just the first order transition 
point. However the corresponding baryon number den¬ 
sity at this point is over 5 times of the normal nuclear 
matter(po = O.I7/m“^). The hadronic EOS is usually 
not reliable at such a high density. Besides, the system 
must be strongly interacting in the transition region, so 
that it can be described neither by an extrapolation of 
the hadronic EOS from the low-density side nor by an ex¬ 
trapolation of the quark EOS from the high-density side 
m- And we should point that it is not a special result in 
our hadronic and quark EOS. Because the LQCD shows 
that that even in relatively large region of chemical po¬ 
tential, QGP is still strongly interacting. So the pressure 
of the quark EOS tends slowly to the free quark gas in 
the P — plane. That means the pressure increasing 
slowly in the P — iib plane, and the intersection between 
quark matter and hadronic matter usually occurs at high 
chemical potential and high number density if we choose 
a modern model to calculate the EOS for quark matter 
rather than a MIT bag or a perturbative model. Now, 
we start to construct EOS for hybrid stars in the point 
of view of the smooth crossover. Erom the Refs. [2Qll22] 



FIG. 2: The pressure as a function of the baryon chemical 
potential. The dashed line is the EOS of hadronic matter, the 
solid line represents the EOS of quark matter. The vertical 
dashed lines show the corresponding baryon number density 
of the hadronic matter. The baryon number density is scaled 
by the baryon number density of normal nuclear matter po = 
O.I7 /?ti“^ 



FIG. 3: The pressure as a function of the baryon number 
density. The shadow region is possible crossover region, the 
solid line is the EOS of quark matter and the dashed line is 
the APR EOS. The x-axis is the baryon number density scaled 
by nuclear number density and the y-axis is the pressure 


we can see that the strategy is not unique, but it should 
meet the requirements of physics. Refer to the Ref. m, 
the so called “3-window-model”, we construct the hybrid 
EOS from the P — pB plane. The pressure as a function 
of the baryon number density is showed in the Eig. 3. 
The shadow region is the possible crossover region that 
is between the 2po and 4po. We adopt the interpolation 
function in the Ref. |20] to make a smooth connection 
between the two EOS in the P — pB plane: 

P = Ph x f- + Pq X /+, 


( 22 ) 
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P(GeV^) 

FIG. 4: The interpolation EOS (green dot line) compared with 
the quark EOS(solid line) and hadronic EOS(dashed line). 


/± = ^ (^1 ± tanh > (23) 

where Ph and Pq are the pressure in the hadronic matter 
and that in the quark matter,The window p — T < p < 
p + r characterizes the crossover region in which both 
hadrons and quarks are strongly interacting, so that nei¬ 
ther pure hadronic EOS nor pure quark EOS are reliable. 
From the thermodynamical relation, P = p^d{elp) jdp^ 
we obtain 

e{p) = eH{p)f-{p) + eQ{p)f+{p) + Ae (24) 

Ae = p r {enip') - eQ{p')) ^dp' (25) 
J P P 

with g{p) = ^ (e^ + e“^) ^ and X = Here eni^q) 
is the energy density obtained from Hadronic EOS and 
quark EOS. Ae is an extra term which guarantees the 
thermodynamic consistency. Correspondingly, if we start 
from the energy density and deduce the pressure (This is 
called e — interpolation in the Ref. EOl), there will be a 
AP term in the expression of the pressure. But we find 
the AP term often leads to a fluctuation in the hybrid 

EOS and makes the sound velocity {vg = become 

larger than light. This is obviously unreasonable, so we 
don’t try this. The final result of this hybrid EOS is 
showed in the Fig. 4, compared with the quark EOS and 
hadronic EOS. From the Fig. 4, the relation between en¬ 
ergy density and pressure, we can evaluate the stiffness 
of the EOS easily. The APR EOS has a well-known stiff 
tail in the large density region. So, in the Fig. 4 we can 
see that when the pressure and energy density are large, 
the slope of the APR EOS becomes much lower than 
that of the DSEs EOS. While, in the low density region, 
the quark EOS is often stiffer than the hadronic EOS. 
The hybrid EOS based on first order phase transition 
equals to the hadronic EOS in low density region, and 



EIG. 5: The pressure as a function of the baryon chemical 
potential. The green dot line shows the pressure of interpo¬ 
lation EOS, the solid line is the pressure of the quark matter 
and the pressure of the hadronic matter is the dashed line. 


in the large density region it equals to the quark EOS. 
So, the hybrid EOS based on first order phase transi¬ 
tion is often softer than both the hadronic EOS and the 
quark EOS. However, if we construct the hybrid EOS 
from smooth crossover, the EOS will be influenced by 
the quark EOS in the low density region, and be influ¬ 
enced by the hadronic EOS in the large density region. 
So, the result will be different. Unfortunately, our hybrid 
EOS is still very soft if we utilize this interpolation strat¬ 
egy, because our quark EOS is soft. This is in accordance 
with the conclusion in the Ref. [20] that the maximum 
mass of hybrid stars can exceed two solar mass only if 
the quark matter has a stiff equation of state and the 
crossover takes place at around three times the normal 
nuclear matter density. However, we find another way to 
construct two solar mass hybrid stars with even a soft 
quark EOS. We apply the same interpolation function to 
the P — pb plane and then calculate energy density from 
thermodynamical relations. The pressure, energy density 
and the EOS are showed in the Fig. 5, Fig. 6 and Fig. 7. 
The chemical potential in the crossover region is between 
1.4 GeV to 3.4 GeV because the first order phase tran¬ 
sition point is at around 2.4 GeV. In order to quantify 
the stiffness of EOSs, the sound velocities are showed in 
the Fig. 8 and Fig. 9. A comparison can be found in 
Ref. [21]. The EOS obtained by the interpolation from 
the P — Pb plane is soft while the EOS obtained by the 
interpolation from the P — pb is much stiffer. From Fig. 
8 we can see that the sound velocity is always smaller 
than 0.4 times of the velocity of light while the sound 
velocity in Fig. 9 is much larger. 

In the P — Pb plane, we can clearly find the differ¬ 
ence between the phase transition of smooth crossover 
and Maxwell construction. Although the pressure tends 
to quark matter in the high chemical potential region 
and tends to the hadronic matter in the low chemical po¬ 
tential region, it is different from both the quark matter 
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FIG. 6: The energy density as a function of the baryon chem¬ 
ical potential. The green dot line shows the energy density 
of the interpolation EOS, the solid line is the energy density 
of the quark matter and the energy density of the hadronic 
matter is the dashed line. 
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FIG. 9: The sound velocity of the hybrid EOS obtained by 
the interpolation from the P — fis plane 



P(GeV)^ 

FIG. 7: The interpolation EOS (green dot line) compared with 
the quark EOS(solid line) and hadronic EOS(dashed line). 
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FIG. 8: The sound velocity of the hybrid EOS obtained by 
the interpolation from the P — ps plane 



FIG. 10: The m-r relation of the hybrid stars based on the 
smooth crossover EOS. The maximum mass is over two solar 
mass 


and the hadronic matter in the crossover region. Finally, 
by integrating the TOV equations with the EOS, we get 
the mass-radius relation of the hybrid stars. It is showed 
in the Fig. 10. The maximum mass is 2.35 times the 
solar mass. There are also other interpolation functions 
to make the interpolation EOS. For example, a polyno¬ 
mial function is used in the Ref. [22]. But we find that 
different interpolation functions cannot make apprecia¬ 
ble difference. If we adopt the polynomial function to 
construct the EOS, the difference of the maximum mass 
is less than 5 percent. This result is apprehensible, be¬ 
cause the interpolation functions should be smooth at the 
boundaries of the interpolating interval. Thus, the val¬ 
ues of different interpolation functions in the crossover 
region do not differ widely. 
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IV. DISCUSSION 

In this paper we introduce our quark EOS with 3- 
flavor quark based on the framework of DSEs to calculate 
the structure of hybrid star. Eor the hadronic phase we 
adopted the APR EOS and we calculate the mass-radius 
relationship of the hybrid stars in the point of view of the 
smooth crossover phase transition from hadronic matter 
to quark matter. The common belief is that the EOS for 
quark phase should be very stiff to construct a two solar 
mass hybrid star. Eor example, in the Ref. [20], the au¬ 
thors performed the interpolation in the P — p plane and 
claimed that no matter what kind of hadronic EOS is 
adopted, the maximum mass of neutron stars can exceed 
2 solar mass on condition that the crossover takes place 
at around three times the normal nuclear matter den¬ 
sity and the quark matter is strongly interacting in the 
crossover region and has stiff EOS. Nonetheless, we find 
that if we start from the P — p plane, the interpolation 
function can generate a stiff EOS and finally construct 


a massive hybrid star compatible with two solar mass, 
even though our quark EOS is relatively soft. 

Eor simplicity, the hyperons are not included because 
the interaction of hyperons with nucleons and among 
themselves are not well determined. The consideration 
of hyperon will soften the EOS for hadronic matter, but 
in our result, the maximum mass of the hybrid stars is 
larger than the maximum mass of both the pure quark 
stars and the pure hadronic stars. So, to what extent 
the influence of hyperons will be is model dependent and 
need further study. 
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